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ABSTRACT 

The traditional method of data presentation for heat exchanger 
surfaces does not permit comparison of individual surface types in 
any simple manner, nor is the published experimental data plentiful 
enough to conduct a meaningful study of the effect of varying 
geometries on surface performance. This data is most commonly 
presented in terms of heat transfer coefficients and friction factors 
referenced to the exposed area as a function of Reynolds number 
based on the minimum free flow area. There are six basic geometric 
parameters for a triangular pitch circular finned tube bank: tube 

diameter; tube transverse and diagonal pitch; fin height, thickness, 
and spacing. 

Baskerville (3) developed a method of surface comparison for 
circular finned tube heat exchangers. Biery (4) developed a method 
for prediction of heat transfer for finned tube banks in crossflow. 
Jameson, et al., (5) developed an empirical relation for correlation 
of pressure drop data across triangular pitch circular finned tube 
banks . 

The correlations for heat transfer and pressure drop for trian- 
gular pitch circular finned tube banks in crossflow are used with 
the comparison method to study the effect of varying geometry on 
surface performance. Optimal fin spacing and height are found for 
specific fin thicknesses, tube diameters, and tube pitch. Performance 
trends for the six possible geometric parameters are established. 

Thesis Supervisor: Professor Warren M. Rohsenow 

Title: Professor of Mechanical Engineering 
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NOMENCLATURE 



Dimensionless factor; defined by equation (22b). 

Minimum free flow dimension between tubes measured 
diagonally between adjacent rows from base to base as 
if fins were not present; see Figure 3, [FT]. 

Minimum free flow dimension between tubes measured 
transversely to the direction of flow from base to base 
as if fins were not present; see Figure 3, [FT]. 

Heat transfer area of the bare base tubular surfaces; 

2 

equals length times heated perimeter [FT ]. 

Exchanger minimum flow area for gas: (minimum area 

can be transverse to flow or diagonal to flow; see 
Figure 3). [FT^]. 

Total fin area of a tube of length L; defined by 

2 

equation (A14), [FT ]. 

Exchanger flow passage frontal area ignoring any 

2 

enhancing surfaces; see Figure 3, [FT ]. 

Ratio of fin area to total transfer area; defined by 
equation (A15). 

Total air heat transfer surface area (same as total 

surface area exposed to gas crossflow), including fin 

2 

surface and edges of fins [FT ]. 



Longitudinal center-to-center tube spacing; see Figure 3, 
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Constant used in the computation of friction factor (f) 
for flow over smooth tube banks; defined by equation 
(32), 


Ca 


Constant required for the computation of radial fin 
efficiency; defined in Figure 5. 




Dimensionless coefficient defined in equation (23). 


c 

p 


Specific heat at constant pressure, [BTU/LBM-^F] . 


s 


Constant used in calculation of friction factor for gas 
flow over smooth tube banks; defined by equation (31). 


D 


Constant required for the computation of radial fin 
efficiency; defined in Figure 5. 




Equivalent diameter defined by equation (20). 


d 

e 


Equivalent diameter of a finned tube; defined by 

equation ( A18) . 


°f 


Fin diameter, [FT]. 


°j 


Equivalent diameter of a triangular pitch finned tube 
bundle; defined by equation (26). 


°h 


Hydraulic diameter; defined by equation (la), [FT]. 




Ratio of hydraulic diameter to the equivalent diameter 
for an equilateral triangular pitch infinite smooth tube 
bank (ETP-I-STB); defined by equation (21). 




Nominal diameter; defined by equation (1b), [FT]. 



J; 




1'. 













8 



D Tube outside diameter, [FTI. 

O 7 L J 

E Factor defined by equation (20). 



ETP-I-STB 

F 



Equilateral triangular pitch infinite smooth tube bank. 

Transverse center-to-center tube spacing; see Figure 3, 
[FT]. 



f 



Friction factor predicted by the Jameson (5) correlation 
and based on total area (A-j.); defined by equation 
(A23) . 



f 



Friction factor based on total area (A-p); defined by 
equation (4a) . 



f 



J 



Friction factor used by Jameson (5); defined by equation 
(29). 



F. 

I n 



Number of fins per inch; defined by equation (46). 

Friction factor based on base area defined by 

equation (4b ) . 



f Friction factor based on total area (A-^) for a smooth 

s T 

surface, defined by equations (31) and (32), 
g^ 32.174 [LBM/LBF] [FT/SEC^]. 

G Mass flux based on minimum free flow area, [A ]; 

c 2 ^ 

defined by equation (2a), [LBM/HR-FT ]. 



G 

G 



f 

N 



Dimensionless geometry factor; defined by equation (A7). 

Mass flux based on free flow area, [Ap]; defined by 
equation (2b), [LBM/HR-FT^] . 
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Heat transfer coefficient based on total area, [A^]; 



defined by equation (5a), [ BTU/HR-FT^-°F ] . 



N 



Heat transfer coefficient based on base area, [A, 1: 

9 n ’ b-* ’ 

defined by equation (5b), [BTU/HR-FT -°F]. 



I 



k 



Projected perimeter of a finned tube; defined by 
equation A17, [FT], 



j 



Colburn J-factor based on total area, [A^]; defined by 
equation (7a). 



■^'b&Y Colburn J-factor based on outside tube diameter [D^]; 

defined by equation (A22). 

Jq Colburn j-factor predicted by Biery's [4] method as 

outlined in Chapter 3. 

j|^ Colburn j-factor based on base area, [A|^]; defined by 

equation ( 7b ) . 



V&L 



Colburn j-factor based on experimental data from Kays & 
London [2] . 



j 



M 



Colburn j-factor based on equivalent diameter, 
defined by equation (A20). 



j 



s 



Colburn j-factor based on total area, [A^], for a smooth 
surface, defined by equation (24). 



k 



Thermal conductivity [BTU/HR-FT-^F] . 



I 



Fin height defined by ^or rectangular fin sheets, 

see Figure 5, [FT ] . 



10 



L 


Flow length along the axis of the tubes (i.e. for the 
fluid inside the tubes). 


m 


Component of fin efficiency ^ ^ f defined by equation 

(10). 


N 


Number of tube rows in a tube bank. 




Parameter used in algorithm of Figure A4 for calculating 
when minimum free flow area occurs in the diagonal 
space. 


*^HT 


Parameter used in algorithmof Figure A4 for calculating 

D. when minimum free flow area occurs in the transverse 
h 

space. 




Parameter used in algorithm of Figure A4 for calculating 
nominal diameter [Oj^]. 


^PR 


Symbol used in Kays and London [2] for PR. 




Symbol used in Kays and London [2] for RE. 


^ST 


Stanton Number [h/G^Cp]. 


Nu 


Nusselt number; defined by equation (6a). 




Nusselt number; defined by equation (6b). 


NTU 


Number of transfer units. 


’^f-corr 


Corrected fin efficiency per Figure 5. 


’'f-L 


Fin temperature effectiveness of a longitudinal fin of 
rectangular profile; defined by (9). 
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^’^corr Correction factor applied to permit the use 

of the simpler [Tanh ml]/ml equation for fin temperature 
effectiveness in lieu of Bessel Function type equations. 
See Figure 4. 

Total surface temperature effectiveness, defined by (8). 

P Pumping Power, [HP]. 



P^ Property constant for fin material and air properties; 

defined by equation (A28). 

2 

AP^ Friction pressure drop, [LBF/FT ], 

Pr Prandtl Number [/^C /K] . 

P 

q Heaf transfer rate, [BTU/HR]. 



RE Reynolds number based on minimum free flow area, [A^] ; 

defined by equation (3a) and labeled Np by ref [2]. 

RE^ Reynolds number based on the equivalent diameter of a 

finned tube handle, mass flux based on 

the minimum free flow area, defined by equation 

(28). 



RE 



N 



r 

e 




Reynolds number based on free flow area, [Ap]; defined 
by equation (3b). 

Radius of enhanced (finned) surface, [FT]. 

Dimensionless geometry factor; defined by equation (A5). 

Hydraulic radius, equals hydraulic diameter divided by 
four for a round tube, [FT]. 
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Outside radius of base tube, [FT]. 



BF 



Ratio of spacing between fins to tube outside diameter. 



Temperature, [^F] . 



T . 



c, I n 



Cold gas inlet temperature, [T . ] 



Ratio of fin thickness to tube outside diameter, [ 8 /D ] 

o 



H, i n 



Hot fluid inlet temperature, [^F]. 



V 



Volume, [FT'^j 



V 



os 



Volume one side; volume of the gas side of the heat 
exchanger on a per tube basis. [FT ] 



V. 



Total volume of heat exchanger on a per tube basis 

[FT^]. 



w 



Mass flow rate, [LBM/HR]. 



X,Y,Z Principle dimensions of heat exchanger, [FT] 



X, 



Ratio of diagonal tube spacing to tube outside diameter. 



X. 



Ratio of longitudinal tube spacing to tube outside 
diameter, equals b divided by D . 



X 



MAX 



The greater of X.^ or X^. 



<\>A 



Ratio of longitudinal spacing to tube outside diameter 
for an equilateral pitch tube bank. 
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(Xj)^ Ratio of transverse tube spacing to tube outside diameter 

for an equilateral triangular pitch tube bank. 

Miscel laneous 



a 


In tube geometry, angle defined by the longitudinal 

centerline and the diagonal centerline; see Figure 4, 
[degrees] . 


/? 


Ratio of total transfer area on one side of the exchanger 
to total volume of the exchanger, labeled as a in reference 

[2], [FtVfT^]. 


8 


Fin thickness, or average fin thickness if tapered, [FT]. 




Viscosi ty , [LBM/HR-FT ] . 


P 


Density, [LBM/FT^] . 



Subscripts 



d 


Diagonal 


e 


Enhanced surface 


ETP 


Equilateral triangular pitch. 


MAX 


The maximum value a variable may assume. 


m 


Heat exchanger fin material 


o 


Outside 


os 


One side 


P 


Predicted 
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REF Reference value for indicated parameter, 

s Smoth surface 

T Total 

1 Tube side; hotter fluid flowing inside smooth tubes. 

2 Finned side, colder fluid flowing over finned 

banks . 



tube 
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I • I ntroduct Ion 
A. Purpose 

The objective of any heat exchanger comparison is to enable 
the designer to select that unique surface which can be termed as 
the most beneficial to his particular application. The comparison 
technique should be logical, easily implemented, and inherently 
accurate. 

Enhanced surfaces, in the form of plate fins and finned 
tubes, have long been used in heat exchanger design to improve 
heat transfer character istics . Many have presented methods of 
comparing the performance of various types of surface designs. 

Soland, et al., (1) developed a method of comparison that 
permits performance comparisons of all the Kays and London (2) 
plate finned and unfinned surfaces. Baskerville (3) developed a 
method of comparison, using theory and principle inherent in 
Soland*s method as a basis, that permits performance comparisons 
of crossflow finned tube surfaces presented in Kays and London. 
Biery (4) developed a method to predict heat transfer coefficients 
for gas flow normal to smooth and finned tubes. Jameson, et al (5) 
developed an empirical relation that correlates pressure drop data 
for circular finned tubes. 

The purpose of this paper is to: 

1. Employ the methods and correlations developed by Biery, 
Jameson, and Baskerville to predict the performance of a 
crossflow circular finned tube surface for any given 



geometry . 
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2. Evaluate the accuracy of the predicted performance by 
comparison to surfaces presented in Kays and London. 

3. Use the performance predictions to evaluate the effect of 
varying geometry on the performance of crossflow 
circular finned tube surfaces. 

4. Develop graphs and tables that will assist the designer 

in selecting the geometry of a crossflow circular finned 
tube surface which will be the most beneficial for his 

appi ication . 

B. Background 

Crossflow finned tubular heat exchangers have found 
applications as gas turbine plant i n tercoolers, aircraft engine 
and electronic coolers, air conditioning units, and various other 
heat exchanger uses. These applications most often involve 
gas-liquid service. [See Figure l]. 

The performance of a given heat exchanger in which the 
ratio of the heat transfer coefficients between the two working 
fluids is greater than about three is markedly improved by the 
employment of enhanced surfaces. The increased surface area 
inherent with the employment of fins yields a net improvement in 

heat transfer despite lower heat transfer coefficients. 

Optimum design is often used synonomously with "lowest 

cost". Cost is frequently defined in terms of amortized acquisition 
cost, operating cost, or impact cost in terms of either weight or 
volume. An example of the latter might be aircraft engine 

component design where weight and volume are valuable 
quantities which can be directly equated to dollar cost. 
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Estimation of dollar cost of a particular heat exchanger 
design is difficult, primarily due to the proprietary nature of the 
required cost estimating relations. The term ’’cost" or "optimum”, 
therefore, as used in this paper, will most generally refer to 
"impact cost" and will serve as a means of relative ranking of 
proposed surfaces specifically by required volume, pumping power 
and heat transfer (NTU) comparison. 
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figure I 





24 



I I , Comparison Method 



Traditional ly , 


data is presented in 


terms of 


heat transfer 


Col burn 


modul US, 


j, and friction 


factor, 


f, based 


on the total 


exposed 


area, A^, 


as a function of 


Reynol 


ds number. 


RE, based on 



minimum free flow area, A^, and a hydraulic diameter, D^, of the 
flow passage. [See Figure 2] 

Baskervi I le*s (3) comparison method converts these j and f 
magnitudes to new quantities and f^^ based on the heat transfer 
area of the bare base tube surface, A^, and a Reynolds number, 
REj^, which is referenced to the open flow passage area, Ap, as 
though fins were not present. The effect of the fins is accounted for 
as an increased heat flux and hence larger h for the bare base 
tube surface area. In order to Incorporate the effect of the fins 
into the metal conductivity must be specified. 

To convert data presented in the format of Kays and London 
(2) to the new basis applicable to the crossflow finned tube case, 
Baskerville (3) derived various ratios from basic heat transfer 
definitions and theory presented by Soland (1). Important assump- 
tions are that the heat transfer resistance of the tube walls is 
negligible and that the controlling heat transfer resistance is 
assumed to be on the finned side of the heat exchanger. Comparison 
assumes that the flow rate, W, the hot fluid inlet temperature, 

T, . and the cold gas inlet temperature, T . , are held 

h , I n c, I n 

constant. Basic definitions are presented in Table I. Ratios for the 
two cases of minimum free flow area which may occur [See Figures 3 
and 4] are presented in Table II. Fin efficiency is calculated per 



Figure 5. 
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FIGURE 2: Sample data plot 



showing trad i t ion a I 



method of data 



Presentation . 




FIG. 92 



flNNED CIROJUJl TUBES 
SURFACE Cf - 

T«W oal«i4« dlwMtar - O.J0 in. 

ria plt«ia - T.34 p«r inch 

Flov hydirMilic - 4r.-0.0154 ft. 

Pin ihiokn4«« (&T«rag«)* - 0.018 in. 

Pr«#>fiow ar««/frontnl •»<|'«0.538 . . 

B«*t trajLsf «r ftran/total toIiumi - or “140 ft. /ft. 

Pin aran/total araa - 0.892 

Rota I Expariaaotal uncartainty for baat tmnafar 
raaalta poaaibly aoaavhat groatar than tha 
BooUaa.1 *5% quoted for tha other aarfacaa 
bacauaa of tha oacaaaitF of aatisating a 
eootaet raaiatanca in tba bi>«atal tnbaa. 



a - Pina alicbtlj taparad. 



26 



LU 

GO 

Z3 



UJ 

cc 

< 

GO 



n 

CM 



< 

II 

2 

O 



iO 

ro 



Q 

2 

O 



UJ 

cc 



n 

<T 



CM 



CL 



O 

iO 

<!■ 



O 

U) 

CL 

CM 



lO 



-Q 



< 

C7 



2 

Q 



II 

2 

D 



CO 

z 

O 



Ll. 

UJ 

Q 



UJ 

CO 

< 



CM 



> 

< 



fO 

CM 



fa 

ro 



(a 

<r 



o 



cc 



(a 

lO 



(a 











U 


0 




o 








CM 


0) 




Q 








o 


Q. 1 — 1 




Z 










CM 




o in 










< 


jc 


-j 




d 


M- 








o 


H 


CO 


0. 




1 ^ 0 


< j: 


Q < 


< ^ 


o 


<1 


iU 


b P- 


CJ 


2 




< 










< 




O 






cr 


















CO II 


1 


II )i 






II 



II 

r 



c_ 

<u 

<u 

E 





(a 


> 






•4-^ 




Q 








o 


> 


O 


0 






(U 


— 


3 


> 


H 


fO 




Z 


C_ 


if) 


< 


XI 


(/) 




> 


la 


a 


X 


:e 



u 






<D 


c_ 




n 

E 

D 

Z 


0 

u 

(0 






ll 


C M- 


(/) 




't. 


X 


c 


^ o 


*0 


0 


^ CJ 


c 


u 




> 


fa 


(D 


'Z 


OJ 


CC 


U- 


X 



Nusselt Number 



QUANTITY KAYS AND LONDON [2] BARE TUBE 



27 




TABLE II: Derived Ratios used to Convert Data Found in the Literature to the New Basis 
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FIGURE 3: The Two Cases of Minimum Free Flow Area Occurence and 
the Associated Geometric Relations of Each Case. 
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CASE I : Minimum Free flow area Transverse to Flow Direction 
Total Volume on one Side: 



V 

os 

L 



dimension into paper along flow length of tubes. 



"A|_b" Volume: 



Aj-b = [X_D - D ] L X, D 
F '■To O'" L c 



Ratio: 



Apb = [X^ - 1] X|_ 

V l<Ix; 

OS T L ' 




CASE I I : Minimum Free Flow Area in the Diagonals 

Total Volume on one Side: 

V = [D ^ X^^ sina- ttD Va] 
os o D o 

2 





"Apb" Volume: 



Apb = 2X^ [Xq - I] 



Ratio: 



Apb = 4[Xp - 1]X|_ 

V X_^sina — ir/4 

os D 




FIGURE 4 : Calculation of [A^bj/V^g for both Case I and II. 
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FIGURE 5: Computing fin efficiency of a Radial fin; Ref (3)- 
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Further, Baskerville (3) shows that for any crossflow heat 
exchanger, where for the same fluid at the same temperature, 
jn and p are constant, the fluid pumping power per unit volume 
on one side is given by: 



V 



os 



W AP, 



V P 
os 



2 

g P^D^\ 
o N os 



D \ 

N os 



(17) 



And for the flow rate, W, also constant: 
NTU A, h,. 4/4 



V 



os 



A, h. , 

b N 

V WC 
os p 






D ^ V 
n os 



( 18 ) 



D 



2 V 



os 



Using equations (17) and (18), it is possible to construct a 
performance parameteter plot of NTU/Vq^ versus P/V^^ from which 
useful performance comparisons between two heat exchangers for four 
different criteria are obtained. The four different comparisons, 
depicted in Figure 6, predict relative performance in the following 
categories : 

Case A: Same heat exchanger shape and volume. 

Case B: Same volume and pumping power. 

Case C: Same pumping power and number of transfer units. 

Case D: Same volume and number of transfer units. 

Baskerville (3) plotted performance curves of all the finned 
surfaces of Kays and London (2). Figure 7 shows the resulting per- 
formance parameter plot presenting only the highest performance 
curve at each nominal diameter (D^^). The study considered sixteen 



NTU 
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Performance Parameter Curves For Two Hypothetical 
Surfaces Showing Points Used in Sample Comparisons. 



FIGURE 6: 



NTU 
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FIGURE 7 

PERFORMANCE PARAMETER CURVE FOR FINNED TUBULAR SURFACES 
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different surfaces having eleven different nominal diameters. Table 
III lists all of the Kays and London surfaces considered by type 
and surface designation in addition to the one surface supplied by 



the Trane Corporation. 




Baskerville (3) also 


offered a procedure which can be used 


for sizing crossflow finned 


tubular heat exchangers. Two designs of 


an intercooler for a 5000 


SHP gas turbine plant were undertaken 



using surfaces 98C and 94. Table IV lists the fluid properties 
assumed for both designs. The principle dimensions using surface 
98C and then surface 94 were calculated to be: 



Surface 98C : 


Surface 94: 


X=6.43 ft 


X=1 .5 ft 


Y=22.11 ft 


Y=10.89 ft 


Z=0.63 ft 


2=1 .76 ft 



3 3 

Total Volume = 105.46 ft Total Volume = 28.68 ft 

Both designs satisfy the thermodynamic and geometric con- 
straints. However, use of surface 94 results in a 73% reduction in 
volume, and is seen from Figure 7 to be notably superior to surface 



98C. 



Surface Numbered Kays and London (2) 
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TABLE IV 



5000 


SHP Gas Turbine Plant Intercooler Fluid Properties 



Propert ies 


Side 1 Side 1 1 


Fluid 


Water Dry Air 


W 


400,000 Ibm/hr 200,667 Ibm/hr 


T. 

in 


60° F 260° F 


"'"out 


82° F 80° F 


P 


1.554 PSI .3060 PSI 
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A performance parameter plot is constructed (Figure 8) for 
surfaces 98C and 94 to quantitatively predict the volume savings. 
This comparison corresponds to Case C. From the design point of the 
reference heat exhchanger (surface 98C) a line of slope 1 Is 
constructed. The intersection of this line with the surface 94 will 
yield the performance parameters at the design point of the new 
heat exchanger. The predicted volume reduction is found from the 
ratio of either the ordinates or abscissas of the two surfaces. 

The performance parameters of the reference heat exchanger 
design (surface 98C) are: (See Figure 8) 

NTU = 401.20 (1/in^) 

V 

OS 

P = 2.93 10^^ (1/in^) 

V 

OS 

The performance parameters that would result using surface 94 
are: (See Figure 8) 

NTU = 1,634.68 (l/in^) 

V 

OS 

P = 1.19 10^^ (1/in^) 

V 

OS 

Since pumping power and number of transfer units are the 
same for both surfaces it is seen that: 




Therefore, a predicted volume reduction of 75% is anticipated 
using the comparison method proposed by Baskervllle. This compares 
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CM 

C 



o 



> 

ZD 




P/V (1/in^) 

OS 



FIGURE 8 



Performance Parameter Plot for the two surfaces used in the 
Design Example. 
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quite favorably with the heat exchanger design calculations made 
for these surfaces and therefore shows that the performance 
parameters for a surface allow comparison of surfaces without going 
through a complete heat exchanger design calculation. Where there 
may well be hundreds of different surfaces to be considered, such 
as the parametric study of geometry effects on crossflow circular 
finned tube heat exchangers conducted in this paper, the value of 
such a comparison technique cannot be overestimated. 
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ill. Heat Transfer Correlations 



Biery (4) presents a method to predict heat transfer 

coefficients for gas flow normal to finned tube banks with 

triangular pitch. A transformat ion from the actual tube bank to an 
equivalent equilateral triangular pitch infinite smooth tube bank 
(ETP-I-STB) is made and the heat transfer coefficient for this 
ETP-I-STB is that desired for the real tube bank. 

The correlation for the desired heat transfer Colburn modulus, 
jj is given as: 

0/7 A 

(19) 

where RE is based upon the hydraulic diameter of Kays and London 
(2) and is given by: 






D. = 4r. = 4A b 

h h c 



(la) 



The minimum flow area, A , often will be in the transverse 

direction to flow. However, for tube banks with a small 

compared to X^, the minimum flow area may well be along the 

diagonal of the triangular pitched bank. Expressions for A^ and 

are developed in Appendix I for both cases. 

Clj is a function of the transverse distance between tubes, 

X^D^, and an equivalent diameter, D^, given by 

= D + E (D. - D ) (20a) 

E o f o 

where for <.04775 inches 

BF o 

E = 1 (20b) 



and for 



.04775 inches— S qi_D —.4 inches 
BF o 

E = 1 ,1357 - 2.842 

Dr O 



(20c) 
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and 



for S_(_D ^.4 inches 

br o 

E = 0 



(20d) 



where SgpD^ is the spacing between fins in inches. Since E is equal 
to one for SgpD^ less than .04775 inches and equal to zero for 
SgpD^ greater than .4 inches, Dp has the characteristic of being 
equal to the fin diameter, D^, for fins with narrow spacing and 
equal to the bare tube diameter, D^, for tubes with either no fins 
or with fins with a wide spacing. 

For an ETP-l-STB the equivalent diameter is simply the 
outside diameter of the tube. Therefore, the ratio of its hydraulic 

diameter, to the equivalent diameter is: 

/D^\ 4[(X^)^ -1] (X^)^ 4[(X^)^ -1] ( V7/2) (X^)^ (21) 

7T 7T 

where (X-j. denotes the transverse tube pitch (ratio of transverse 

tube spacing to tube outside diameter) for an ETP-I-STB. To 
transform from the actual tube bank to an equivalent ETP-l-STB, 
the ratio of the hydraulic diameter to the equivalent diameter 
[Dj^/Dp] is equated to the ratio of the hydraulic diameter to the 

equivalent diameter of the ETP-I-STB [(D^/D^)^ ]. then 

solved to yield: 

(X^).= 1 + (4A + 1)^ (22a) 

T A 7 ; 



where 



of 




4 (yT/2) 



Biery found that with 
only (^j)^ 



(22b) 

this transformation , was a function 

property of bringing all of the various 
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conf igurations onto one curve. Correlating the data given by Kays 

and London (2), Jameson (5), Fraas and Ozisik (6), Mirkovic (7), 

and Grimeson (8) with empirical curve fitting resulted in: 

= .2818(X-j.)^ - .1282 - 0.08263 [l.2 - (23) 

L(X^)^ - .8924J 

Biery’s transforma t ion and correlation for heat transfer 
prediction of circular finned tubes is used later in this paper to 
study the geometry effects on the sizing of crossflow circular finned 
tube heat exchangers. It should be noted that the data from which 
this correlation is derived is based on a limited range of tube 

spacing, tube diameters, fin spacing, fin thickness, and fin height. 
Application of this correlation beyond this limited range should be 
recognized as an extrapolation of the empirical data. 

A study was made to compare Biery*s method for prediction of 
heat transfer coefficients with other correlations and the data 
presented in Kays and London (2) for circular finned tubes. The 
details of this study are presented in Appendix II. In general it 
was found that Biery’s method resulted in heat transfer coefficients 
that agreed with Kays and London data within 20% of the mean 
value with a mean 7% below the Kays and London data and a 
standard deviation of 8%. 

Zhukauskas (9) developed a correlation for heat transfer for 
turbulent flow over staggered tube banks with the appropriate 
nominal diameter, D.,, which is also used later in this paper to 
study geometry effects on the sizing of crossflow circular finned 
tube heat exchangers. This correlation for smooth tube bank Colburn 
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Modulus, is given as: 



• - 7c; iKi \-0.396 1^ j/6 (D /D.,) 

= .35 (Npj^) (Xj/X^) o' W 



.400 



(24) 



To justify the use of the above relation it is noted that 
Baskerville (3) used this relation to predict the heat transfer for 
crossflow over bare tube banks presented in Kays and London (2). 
The results showed a maximum deviation between the experimental 
value of from Kays and London (2) and the calculated value of 
from equation (24) of 7.5%. These results are shown in Figure 9. 



FIGURE 9 

Deviation of the Zhukauskas (9) correlation for forced- 
convection, turbulent flow over staggered tube banks 
in predicting experimentally obtained value of Colburn 
modulus [j ] as presented in Kays and London (2). 
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I V . Pressure Drop Correlations 

In general the friction factor Is a function of the Reynolds 
number : 

f = C^RE (25) 

where and are constants, 

Jameson (5) developed the following empirical relation for an 
’’equivalent diameter”, Dj, of a triangular pitch finned tube bundle 
which correlated his data for pressure drop across tube banks in 
which the minimum free flow area is in the space between tubes 
transverse to the air flow. 



D. 

J 




(26) 



where d is defined by equation (A18). 



Using this expression for Dj, equation (25) becomes: 
f = 1 .532 RE 

J 



(27) 



where REp. is the Reynolds number based on D. and G . 

U I J c 



RE 



G D, 
c J 



(28) 



D 



The friction factor is defined by Jameson (5) for N tube rows 
as follows: 

O 2f,G 
= J c 



P 9. 



( 29 ) 
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This definition differs from the Kays and London (2) definition 
which is given by: 






ffG ^N1 
c 




} 

1 > 
|H 

1 




'2fG 1 

c 




[\°ol 


'■ 9oJ 








L P9o j 




L J 



(30) 



Equation 


(30) 


is the same as 


equation 


(4a) , 


except 


that 


equation (4a) 


appi 


ies to a single 


row of 


tubes. 


where It 


is 


understood that 


the 


pressure drop is 


direct 1 y 


proportional to 


the 


number of tube 


rows . 













To compare the friction factor used in the Jameson correlation, 
fj, with the friction factor given in the data by Kays and London, 

the geometry factor, (X. D /D^^), must be included and the Reynolds 

number converted to a common diameter. This is done in Appendix 

ill. In general, it was found that the Jameson correlation predicted 
the Kays and London data to within approximately 20% from the 
mean value with a mean of +7,8% and a standard deviation of 7.3%, 
As noted by Jameson, this correlation does not accurately predict 
the friction factor for finned tubes in which the minimum free flow 
area occurs in the diagonal of the triangular pitch tube bank. Not 
withstanding this difficulty, Gianolio and Cute (11) in their set of 

experiments for pressure-drop measurements across equilateral trian- 
gular pitch finned tubes banks in induced draft, noted the Jameson 
correlation compared the most favorably with their results showing 
a maximum deviation of only 15% compared to their data. In a 
study of the geometry effects on triangular pitch finned tube banks, 
the effect of tube pitch can be demonstrated with smooth tubes. The 
other geometry effects can then accurately be demonstrated by 



48 



selecting a tube pitch in which the minimum free flow area does not 
occur in the diagonal. Therefore, the Jameson correlation is used 
later to study geometry effects by predicting performance parameters 
for tube banks in which the minimum free flow area occurs in the 
spaces transverse to the air flow. 

To predict the friction factor for smooth tubes, (^g)> Basker- 
ville (3) derived the following empirical correlations based on the 
Kays and London (2) data for both cases of minimum free flow 
area. 



For Case I where the minimum free flow area occurs in the 



spaces transverse to flow: 



f. = .7184 [X^] , 



(31) 



where: 



-0.4277 



= .828 - 16.596 RE^ 

For Case II where the minimum free flow area occurs in the 



diagonal spaces: 



(32) 



where: 






- 0.7284 



To justify the use of the above relations, it is noted that 
Baskerville used these relations to predict the actual experimental 
friction factor data for flow over bare tube banks as presented in 
Kays and London (2). The results are shown in Figure 10. The 
maximum deviation between the Kays and London (2) experimental 
value of fj and the calculated value from either (31) or (32) is 
approximately 3.0% from the mean value. 
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V . Comparison of Performance Parameter Prediction Program . 

Appendix IV presents an algorithm and calculator program 
that employs the comparison method introduced in Chapter II and 
the heat transfer and pressure drop correlations introduced in 
Chapters I I I and IV to predict the performance parameters of a 
surface. This program was used to predict the performance 

parameters for the triangular pitch circular finned surfaces in Kays 
and London (2). These performance parameters were then compared 
to the performance parameters calculated by Baskerville (3) from 
the actual heat transfer and pressure drop data given in Kays and 
London (2) for these surfaces. The predicted performance parameters 
and the performance parameters calculated by Baskerville are shown 
graphically in Figure II, which offers a visual comparison between 
the predicted performance parameters and the performance para- 
meters generated using the Kays and London (2) experimental data. 
Each surface was compared on the basis of same pumping power and 
number of transfer units (Case C). With this comparison the ratio 

(V ) /(V ) can be found. The percent deviation in (V ) from 

os p ' os e os p 

(V ) can then be found from this ratio and calculated according 
os e 

to equation (33) . 

(V ) - (V J (33) 

Percent Deviation = — - — x 100% 

os e 

The average deviation of the predicted volume one side 
from the calculated (from experimental data) volume one side 
is given by Table V. If surfaces 93, 94, and 98D are discounted 

due to the difficulties already cited for these surfaces, the average 
deviation for the remaining surfaces is +11.7% with a standard 
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Figure 11. Comparison of predicted performance parameters 
performance parameters calculated by Baskerville 
from Kays and London (2) experimental data. 
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Figure 1 1 . (Cont) 
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TABLE V 

Deviation of performance parameter program of Figure A5 for volume 
one side in predicting experimentally obtained value as presented 
in Kays and London (2) and calculated by Baskerville (3). 

Surface Average Deviation 



92 


+23% 


93 


+86% 


94 


+51% 


96 


+24% 


97A 


+30% 


97B 


+20% 


98A 


+5% 


98B 


-6% 


98C 


-11% 


98D 


+25% 


98E 


-2% 


99A 


+24% 


99B 


+ 10% 
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deviation of 13*8%* This degree of accuracy should provide the 
correct indications in the parametric study of geometry effects on 
triangular pitch circular finned tubes presented in the next 



chapter* 
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V I • Geometry Effects on Triangular Pitch Circular Finned Tube Heat 
Exchangers in Crossflow. 

A . General 

A parametric study was conducted using the program of Figure 
A5 to determine the effect of varying tube pitch, diameter, and fin 
geometry on exchanger effectiveness. Copper fin material and dry 
air at 90 F on the enhanced surface side was assumed. Each 

surface was compared to a reference surface for the same pumping 
power and number of transfer units. This corresponds to Case C 

discussed in Chapter li. Therefore, the relative figure of merit for 
the surface in question is the ratio of the volume one side of the 

reference surface to the volume one side of the surface in question, 
(V . With the volume of the reference surface in the 

numerator, this relative figure of merit can be thought of as the 
number of times larger (or smaller) the heat exchanger composed of 
the reference surface would have to be in order to have the same 
performance as the surface to which it is being compared. The 
reference surface is taken to be a bare tube surface with a one 
inch tube diameter and an equilateral triangular tube pitch of two. 
The reference surface is completely an arbitrary choice. 

B . Effects of Tube Diameter and Tube Pitch 

The performance parameters produced by the program of Figure 

A5 for smooth tubes (no enhancements) were used to study the effect 

of varying tube diameter and tube pitch. Table VI lists the values 

of (V ) /V which were generated for various values of tube 

os REF os 

diameter and tube pitch. These values are plotted on a log-log 



scale in Figure 12. 
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TABLE VI: (V )occ/^ ^or various values of tube diameter and 

os KLr os 

tube pitch for smooth tubes. The reference surface is a 
smooth tube with a one inch tube diameter and equila- 
teral triangular pitch of two. 

A. (V varied tube diameter and several tube pitches. 

os Ktr OS 



X^=1.1, X^=2.5, X^=3.5, 



D ( 1 nches) 
o 


><ETP=2 


Xetp=^ 


X 

o 

II 


X^=1.5 


X^=2.0 


.1 


28.2 


4.65 


316 


81 .7- 


31 .0* 


.3 


5.74 


.929 


- 


- 


- 


.5 


2.74 


.440 


32.4 


7.85- 


2.99* 


.75 


1 .52 


.243 


- 


- 


- 


1 .0 


1 .0 


.159 


12.02 


2.86- 


1 .09* 


2.0 


.366 


.0575 


- 


- 


- 


2.5 


- 


- 


- 


.755- 


.288 


4.0 


.134 


.0208 


- 


- 


- 


5.0 


- 


- 


- 


.275- 


.105 


10.0 


.0354 


.00545 


- 


- 


- 


B. Variation in equilateral triangular tube 


p i tch . 




X 




D =1" 


D = 


.25 




ETP 




o 


o 






1.1 




14.5 


106 




1 .3 




4.83 


35.6 




1 .5 




2.61 


19.3 




2.0 




1 .0 




7.48 




3.0 




.329 




2.49 




5.0 




.0918 




.705 




10.0 




.0172 




.135 




* Indicates 


minimum free flow area occurs 


in the diagonal 



(2X|^< X^+1). 
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TABLE VI. (Cont) 

C . Variation in diagonal tube pitch. 





D =1", D =2.5", 

o’ o ’ 


II 

o 

a 


o 

o 

II 


II 

o 

Q 


Xd 


X 

H 

II 

1 


CM 

II 

1 — 
X 


Xt=3 


II 

1 — 
X 


X 

H 

II 

cn 


1 .12 


12.0* 


91 .2* 


D.N.E. 


D.N.E. 


D.N.E. 


1 .3 


6.21* 


46.7* 


D.N.E. 


D.N.E. 


D.N.E. 


1 .49 


4.02* 


30.1* 


D.N.E. 


D.N.E. 


D.N.E. 


1 .51 


2.08 


15.4 


D.N.E. 


D.N.E. 


D.N.E. 


1 .59 


- 


- 


2.16- 


D.N.E. 


D.N.E. 


1 .75 


1 .36 


10.2 


1 .64- 


D.N.E. 


D.N.E. 


1 .99 


- 


- 


1 .24- 


D.N.E. 


D.N.E. 


2.0 


1 .0 


7.48 


- 


D.N.E. 


D.N.E. 


2.07 


- 


- 


- 


1 .00* 


D.N.E. 


2.25 


- 


- 


- 


.789* 


D.N.E. 


2.49 


- 


- 


- 


.643* 


D.N.E. 


2.55 


- 


- 


- 


- 


.607* 


2.75 


- 


- 


- 


- 


.473* 


2.99 


- 


- 


- 


- 


.401 * 


3.0 


.466 


3.52 


- 


- 


- 


5.0 


.207 


1 .59 


- 


- 


- 


10.0 


.0746 


.585 


- 


- 


- 



^Indicates minimum free flow area occurs in the diagonal 

(2Xp<X.^ + 1) 

D.N.E. = Does Not Exist; indicates geometry not possible for the 
indicated values of tube pitch. 
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TABLE VI . (Cont) 



D. Variation 


in transverse 


tube pitch. 




D =1", X^.,=2 
o ’ D 


D^=1", X^=3 


1 .1 


2.34 


1.11 


1 .5 


1 .36 


.653 


2.0 


1 .0 


.466 


2.5 


.883 


.377 


2.99 


.914 


.329 


3.0 


- 


- 


3.01 


1 .22* 


- 


3.5 


1 .09* 


.305 


3.87 


1.11* 


- 


4.99 


- 


.364 


5.01 


- 


.398* 


5.5 


- 


.385* 


5.91 


— 


.417* 



* Indicates minimum free flow area occurs in the diagonal 

(2Xj^<X^ + 1) 



os'REF' os 
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19A (V ) /V versus D for several values of tube pitch 
^ os^REF^ os o 

for a bare tube surface. 



F i g u re 



os'REF' os 
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Figure 12B. 



^^os^REF^^os 
diameter for 



versus values of tube 

bare tube surface. 



a 



os'REF'' os 
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Figure 12C. (V )„_^/V versus for a bare tube surface for 

^ os REF' os D 

several values of tube diameter and transverse tube 



pitch. 



os'REF' os 
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X 



T 



Figure 12D. ( Vos)r£f/'^qs ^ surface for two 

values of diagonal tube pitch. 
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Table VI and Figure 12 show the tube diameter and the tube 
pitch for Case I, minimum free flow area in the transverse space, 
to be approximately proportional to exchanger volume as follows: 



V “ °o' 




(34) 




1 .5 <X.|.< 3.0 


(35) 


2 1 5 


1 .5 <Xp S3.0 


(36) 


^OS “ ^ETP 


1.5<X^.^p<3.0 


(37) 



This leads to the following relation for Case I, minimum free 
flow area in the transverse space (X^ 2X|^-1 ) 



^OS ci 


"><T 


.85 


"^D 


2.15 


" D 

o 


( 38 ) 

, 1.5SX<3.0 


^^OS^REF 


^^T^REF 




f^REF 









Case II, minimum free flow area in the diagonal (20^5 X^+1), 



is not as well defined. Again Table VI and Figure 12 show for Case 
II that the tube diameter and the diagonal tube pitch is 
approximately proportional to exchanger volume as follows: 

1 .45 



V « D 
OS 0 



V « X 
OS D 



V « X 
OS D 



3.85 



2.5 



1 .1 c X|^< 1 .5 



1 .5 <X|^ <=3.0 



(39) 

(40) 

(41) 



However, Figure 12 shows that exchanger volume even on a 
log-log scale is a non-linear function of the transverse pitch. 
Therefore, variation in volume versus transverse tube pitch on a 
log-log scale does not lend to a linear approximation even for a 
very limited range. However, since the change in volume is small 
over the possible range of transverse pitch to a first approximation 
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volume can be estimated to be independent of transverse tube pitch 
for Case II. This leads to the following approximations for Case II: 

- (42) 



V 




3.85 


“ D 


OS ^ 


D 




o 


^^os^REF ~ 


^^D^REF 




f^REF 








— ^ 



1 . 1 *= X. 



1 .5 



V 



os 



os'REF 





■Xd 1 


2.5 


' D 

o 




^^D^REF 




"^^REF 



-1 1 .45 



, 1.5<Xp-=;3.0 



(43) 



XtV4 4 



2 “ 

+ 1)2<2X^<X^ + 1 



It should be recalled that Case II occurs when 2Xj^ < X,^+1 , 
but that 2X|^ c X^ is a geometric impossibility. As the latter limit is 
approached the tube bank approaches an in-line tube bank. Since 
X|^ must be at least j tube diameter, this further restricts the 
range for Case I I . Since: 

(44) 

Case II is restricted to the following range: 

(45) 

This range of possible values is less than a half tube diameter. 
Due to this limited range of geometries in which Case II may occur, 
the minimum free flow area usually occurs in the transverse space 
for most triangular pitch tube banks. 

Finally, a note should be made concerning the discontinuity 
apparent in Figures 12C and 12D for the transition between Case I 
and Case I I . The curves for each Case are produced by a separate 
set of correlations for heat transfer and pressure drop. Therefore, 
unless both correlations are exactly correct, a perfect match would 
not be expected. The amount of mismatch is in a sense a measure 
of the total error of the four correlations taken as a group at the 
point of transition between Case I and Case II. 
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C . Effects of Fin Geometry 

The performance parameters produced by the program of Figure 

A5 for an enhanced surface with circular fins and equilateral 

triangular pitch were used to study the effects of varying fin 

geometry. Table VII lists the values of (V generated for 

os REF os ^ 

various values of fin thickness, fin height, and fin spacing for two 

tube diameters (.4 inch and 1.0 inch) and two values of equilateral 

triangular pitch (X^^p = 2 and X^^p = 3). These values are plotted 

in Figure 13. Again, (V )ot-r- 'S taken to be a smooth tube with a 

one inch diameter and an equilateral triangular pitch of two. Fin 

spacing is given in terms of fin thickness greater the 

value of (V )^_-/V , the smaller is the V of the surface in 

os REF os’ os 

question for a given number of transfer units and pumping power. 

Figures 13A, 13C, 13E, and 13G show the effect of varying fin 

spacing with all other parameters held constant for several values 

of fin thicknesses. The fin height is taken to be the maximum 

possible value, optional value of fin 

spacing, i.e., a maximum value of ^^o 5 ^REp/^os therefore, a 

minimum value of V for the surface in question, occurs for each 

os 

case. This optional value of fin spacing ranges from .8 to 1.2 fin 
thicknesses. The higher values of optional fin spacing in terms of 
fin thicknesses corresponds to the thinner fins. 

Figures 13A, 13C, 13E, and 13G show for a given fin spacing 
in terms of fin thicknesses the surface performance is better for 
thin fins than for thick fins. This is not necessarily true if fin 
spacing is not expressed in terms of fin thicknesses. Figure 131 
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plots versus fin thickness (t^) for two constant values 

of fin spacing , 

Dr;. Again optimal values are found but here 

thicker fins are seen to produce a surface with better performance 

if the fin thickness is less than the optimal fin thickness. Figure 

13J plots ^^os^REF'^^os ^\n values of fin thickness. 

Again it can be seen that for low values of F. the thicker fin 

I n 

produces the better performing surface. 

The information plotted in Figures 13A, 13C, 13E, and 13G can 

be displayed in the form of Figure 13J, since the number of fins 
per inch is given by: 

^ 1 1 (46) 

(SBF+tf)Do (Sg|-/t^+1 )t^D^ 

There may also be other possible forms of displaying the same 

information . 

Figures 13B, 13D, 13F, and 13H show the effect of varying fin 

height from the maximum possible value for values of fin spacing 
which were found to be optimal in Figures 13A, 13C, 13E, and 13G, 

respectively. An optimal value of fin height occurs ranging from 
40% to 90% of the maximum possible fin height. For thinner fins the 
opti mal value of fin height occurs at lower fin heights and is more 
pronounced for thinner fins. 

For some opti mal values displayed in Figure 13, the optimal 
regions are quite flat. This suggests that in the actual design of a 
surface great savings of fin material might be made by only 

slightly compromising surface performances. 
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TABLE VII. ^^os^REF^^os various values of fin thickness, fin 

height, and fin spacing for two tube diameters (.4 inch 
and 1.0 inch) and two values of equilateral triangular 
pitch (X^^p = 2 and = 3). (V^^)p^p is a bare 

tube with a one inch tube diameter and equilateral 
triangular pitch of two. 

A. (V^^) p^^p/V^^ versus fin spacing for four fin thicknesses 
( t^ = .005, .01, .02, .04) for an enhanced tube surface 



with 


XeTP 


= 2, and 


= (H^: 


> 

X 

II 

Cn 




(V for: 

os REF os 






Sep/'f 


t^=.005 


t^=.01 


t^=.02 


t^=.04 


.3 


22.7 


17.3 


13.2 


9.60 


.4 


24.6 


18.5 


14.1 


10.19 


.6 


26.3 


19.9 


15.0 


10.72 


.8 


26.8 


20.5 


15.3 


10.84= 


.9 


27.0 


20.6 


15.3* 


10.80 


1 .0 


27.2* 


20.7* 


15.2 


10.72 


1 .25 


27.0 


20.5 


14.9 


10.38 


1 .5 


26.9 


20.0 


14.6 


10.07 


2.0 


26.0 


19.0 


13.7 


9.29 


2.5 


24.2 


18.0 


12.6 


8.67 


3.0 


22.6 


16.8 


1 1 .8 


8.11 


3.5 


21 .2 


15.9 


11 .2 


7.60 



Indicates Maximum Value. 
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TABLE VII. (Cont) 



B. ^ '^os ^REF^'^^os '^®'"sus fin height for four fin thicknesses 
(t^= .005, .01, .02, .04) for an enhanced tube surface with 

= 1” and ^£jp = 2 and with fin spacing at the opti mal 
value for H, - = .5 



*'^os’rEF^''os 



for : 



^f 


Wr'-° 


Sgp/.f-I.O 




^BF'^V 


t =.005 


o 

II 


CSJ 

o 

II 

4 — » 


o 

It 


.50 


27.2 


20.7 


15.2 


10.84 


.45 


28.7 


21 .8 


15.8 


11.10* 


.40 


30.4 


22.4 


16.2 


11.08 


.35 


31 .7 


22.9* 


16.3* 


10.76 


.30 


32.7 


22.9 


15.8 


- 


.25 


32.9* 


- 


14.8 


9.03 


.20 


32.0 


19.9 


- 


- 


.15 


29.1 


- 


10.5 


- 


.10 


23.2 


12.5 


- 


4.23 


.05 


13.3 


— 


- 


- 



Indicates maximum value. 



71 



TABLE VI I . (Cont) 



• (V )dcc/^ versus 

os REF os 


f i n 


spacing for 


three fin thicknesses 


o 

o 

o 

II 

4-» 


for 


an enhanced 


tube surface with D 

o 


= = 2, and 




- <»f>MAX - 


.5. 


(V ) /V for: 

os REF os 









'BfV 


= 0.1 


= 0.2 


= 0.4 


.4 


35.9 


27.0 


20.1 


.8 


41 .0 


30.6 


22.4 


.95 


- 


- 


22.4* 


1 .0 


41 .7 


30.9* 


- 


1.15 


41 .8* 


- 


- 


1 .2 


41 .8 


30.9 


22.2 


1 .6 


41 .1 


30.1 


21 .2 


2.0 


40.0 


29.0 


20.2 


3.0 


36.6 


26.0 


17.6 


4.5 


32.0 


22.2 


15.2 



^Indicates maximum value 
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TABLE VII. (Cont) 



D. (V^^) p^p/V^^ versus fin height for three fin thicknesses 

(t, = .01, .02, .04) an enhanced tube surface with D = 

t ' ' o 





.4" and X___ 

t \ r 


= 2 and with fin 


spacing at 




value for = 








^^os^REF^^os 


for: 






SBptf=1.15 


Sep/tf^i .05 


Sap/.,-.95 




t^=.01 


t -.02 


o 

II 


.50 


41 .8 


30.9 


22.4 


.45 


45.1 


33.1 


23.7 


.40 


48.4 


35.1 


24.7 


.35 


51 .4 


36.6 


25.1* 


.30 


53.8 


37.4* 


24.8 


.25 


55.1* 


36.9 


23.6 


.20 


54.3 


34.7 


21 .2 


.15 


49.8 


30.1 


17.7 


.10 


39.7 


22.7 


13.3 



^Indicates maximum value. 
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TABLE VII. (Cent) 



E. 


^^os^REF^^os 


versus 


fin spacing 


for four 


fin thicknesses 




(t^ = .005, 


o 

o 


, .04) 


for 


an enhanced tube surface 




with D = 1", 
o ’ 


^ETP = 


3 and 






1 .0. 




‘''os'ref''''. 


for : 
os 










'Bp/'f 


t^ = .005 




= .01 




t^ = .02 


t. = .04 

T 


.40 


14.75 




1 1 .71 




9.11 


6.92 


.60 


15.66 




12.51 




9.74 


7.32 


.775 


- 




- 




- 


7.38* 


.80 


15.88 




12.75 




9.90 


- 


.825 


15.88* 




12.75* 




9.90* 


- 


ID 

00 


15.88 




12.75 




- 


- 


1 .0 


15.77 




12.70 




9.84 


7.26 


1 .2 


15.48 




12.50 




9.67 


7.07 


1 .6 


14.70 




11.92 




9.17 


6.77 


2.0 


13.83 




11.26 




8.63 


6.44 


2.5 


12.80 




10.46 




7.99 


6.02 


3.5 


11 .03 




9.06 




7.06 


5.30 


4.5 


9.63 




7.96 




6.30 


4.75 



Indicates maximum values. 









s 
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TABLE VI . (Cont) 



'f 


F. (V )Dirir/^ versus 

os REF os 

(t^ = .005, .01, .02, 

with = 1", X^.pp 

optimal value for = 

(V )ocir/^ 
os REF os 

Sgp/t^=.825 Sgp/t^=.825 

t^=.005 t^=.01 


fin height for four fin thicknesses 
.04) for an enhanced tube surface 
= 3 and with fin spacing at the 

S3p/t^=.825 S3p/t^=.775 

t^=.02 t^=.04 


.0 


15.88 


12.75 


9.90 


7.38 


.9 


17.46 


- 


- 


- 


.85 


- 


- 


- 


7.86 


o 

00 


19.09 


14.99 


11 .21 


7.93 


.775 


- 


- 


- 


7.94* 


.75 


- 


- 


- 


7.94 


.70 


20.77 


- 


- 


7.91 


.65 


- 


- 


11.66 


- 


.625 


- 


- 


11 .67^^ 


- 


.60 


22.39 


16.74 


11 .64 


7.51 


.55 


- 


16.96 


11 .48 


- 


.50 


23.75 


17.03- 


- 


- 


.45 


24.21 


16.89 


- 


- 


.42 


24.36 


- 


- 


- 


.40 


24.41* 


16.48 


10.04 


5.70 


.35 


24.24 


- 


- 


- 


.30 


23.52 


- 


- 


- 


o 

CM 


19.57 


9.69 


5.50 


2.85 



*lndicates maximum value. 
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TABLE VI I . (Cont) 



G. (V_)occr/V versus fin spacing for three fin thicknesses 
os Ktr os 23 

(t^ = .01, .02, .04) for an enhanced tube surface with 






(V ) /V for: 
os REF os 





t^=0.1 


t^=.02 


t^=.04 


.4 


20.87 


15.87 


1 1 .97 


.8 


23.67 


17.97 


13.39 


.95 


- 


- 


13.47* 


1 .0 


23.98 


18.17 


13.46 


1 .07 


- 


18.18* 


- 


1.1 


24.02* 


- 


- 


1 .2 


23.98 


18.14 


13.35 


1 .6 


23.51 


17.71 


12.88 


2.0 


22.77 


17.07 


12.28 


3.0 


20.71 


15.35 


10.78 


4.5 


17.99 


13.14 


9.43 



*lndicates maximum value. 
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TABLE VI I . (Cont) 



H* ^^os^REF'^^os height for three fin thicknesses 



(t. 



.01, .02, .04) for an enhanced tube surface with 



Dq - ~ 3, and with fin spacing at the optimal 

value for H, - - 1.0. 



*'^os'rEF^'^os 



for: 





t^=.01 


SBp/t,-1.07 
t -.02 


SBF/tf=.95 

t =.o* 


1 .0 


24.02 


18.18 


13.47 


.9 


26.30 


19.77 


14.50 


.8 


28.74 


21 .31 


15.37 


.7 


31 .23 


22.71 


15.98 


.6 


33.66 


23.77 


16.13* 


.5 


35.68 


24.13* 


15.57 


.4 


36.57* 


23.22 


14.08 . 


.3 


34.87 


20.33 


11.54 


.2 


28.41 


15.04 


8.18 



Indicates maximum value 
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TABLE VII. (Cont) 



I. - versus fin thickness for two values of fin 

os RcF os 

spacing (Sgp = .04, .10) for an enhanced tube surface with 



= 1", 


) Xpyp 2, and 




os 


REF^'^os 








^BF = 


.005 


14.13 


7.96 


.01 


14.74* 


9.15 


.02 


13.42 


9.38* 


.04 


10.67 


8.67 


.06 


8.80 


7.89 


.08 


7.50 


7.21 


.10 


6.56 


6.65 


.20 


4.13 


4.88 


.40 


2.52 


3.39 
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TABLE VII. (Cont) 



(V )ppp/V 
OS REF os 



versus fins per inch for two values of fin 



thicknesses (tp = ,01, .02) for an enhanced tube surface 



with = 1", 



= 2, and Hj = 



= .5. 



(V )pcp/V for: 
os REF os 



^in 


t^=.01 


t^=.02 


2.4 


4.25 


4.98 


4.55 


- 


6.71 


6.25 


7.11 


- 


8.33 


- 


9.79 


11 .8 


10.89 


- 


13.3 


- 


12.30 


16.7 


13.24 


13.65 


20.0 


14.79 


14.62 


25.0 


16.79 


15.21 


28.6 


17.99 


15.31* 


33.3 


18.97 


14.62 


40.0 


19.95 


12.02 


45.5 


20.65 


- 


47.6 


20.90* 


- 


50.0 


20.65 


- 


66.7 


19.28 


- 


80.0 


15.85 


- 



* Indicates maximum value. 



os'REF' os 
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Figure 13A. (V )pcp/'^ versus fin spacing for four fin thicknesses 

OS t. r OS 

(t^ = .005, .01, .02, .04) for an enhanced tube surface 



with D - 1 ' 



^ETP " *^f 



os'REF' os 
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Figure 13B. (V )pcp/V versus fin height for four fin thicknesses 
(t = .005, .01, .02, .04) for an enhanced tube surface 

with D = 1” and = 2 and with fin spacing at the 

optimal value found in Figure 13A for 



os'REF' os 
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Figure 13C. (V )/occ/^ versus fin spacing for three fin thick- 

OS Ktr OS 

nesses (t^ = .01, .02, .04) for an enhanced tube 

surface with = .4", = 2, and H^=(H^)j^^^=.5 



os'REF' os 
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Figure 13D. (V versus 

^ os REF os 

(t^ = .01, .02, .04) 

D = .4" and X,--- 

O c 1 r 

Optimal value found 



fin height for three fin thicknesses 
for an enhanced tube surface with 
= 2 and with fin spacing at the 

in Figure 13C for 



os'REF' os 
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Figure: 



13E. 



(Vos)pEF'^^os spacing 

nesses (t^ = .005, ,01, .02, .04) 

surface with D = 1", X__r, = 3, H 
o ETP 



for four fin thick- 
for an enhanced tube 

f " ^^f^MAX = 1 .0 



OS REF' os 
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Figure 13F. (V )or-ir/'^ versus fin height for four fin thicknesses 
OS REF os 

(t^ = .005, .01, .02, .04) for an enhanced tube surface 

with D = 1", = 3 and with fin spacing at the 

o ETH 

optimal value found in Figure 13E for 



os'REF' os 
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Figure 13G. (V versus fin spacing for three fin thick- 

os REF os 

nesses (t^ = .01, .02, .04) for an enhanced tube 

surface with = .4", ^£-pp = .0 



os'REF' os 
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Figure 13H. 



(V )„^^/\/ versus fin height for three fin thick- 
os REF os 

nesses (t^ = .01, .02, .04) for an enhanced tube 

surface with = .4", ~ 3, and with fin spacing 

at the optimal value found in Figure 13G for = 



<«f>MAX 



1 .0 



os'REF' os 
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Figure 131. 



^^os^REF 




versus fin 


thickness 


for 


two 


values of fin 


spacing 


^^BF 


o 

o 

ti 


for an 


enhanced 


tube surface 


with D 

o 


= 1", 


><Etp 


and = 


(H^) 


MAX 


= .5 



.5 



OS REF' os 
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F. 

I n 

Figure 13J. ( versus fins per inch for two values of fin 

thicknesses (t^ = ,01, .02,) for an enhanced tube 

surface with D = 1", = 2, and H, = = .5 

o ETP f f MAX 
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Effect of tube diameter and pitch for enhanced tube surfaces 

can also be seen from Figure 13 from the two values of tube 

diameter and the two values of tube pitch for which (V )_--/V 

^ os REF os 

was plotted for the various fin geometries. The effect of tube 
diameter and tube pitch on exchanger volume for an enhanced 

surface is less than a bare tube surface. If the optimal points in 
Figure 13 are compared for the two tube diameters and two tube 
pitches, the is approximately doubled for a 2.5 increase in 

tube diameter and increased 1.5 times for a 1.5 increase in 

equilateral triangular tube pitch. For a bare tube surface larger 
changes in would be predicted for the same change in tube 

diameter and equilateral triangular tube pitch. 

D . Summary, Conclusions, and Recommendations for Future Studies. 

Correlations for heat transfer and pressure drop for trian- 
gular pitch circular finned tube banks in crossflow developed by 

Biery (4) and Jameson (5) are used with a method for comparison of 

finned tube surfaces developed by Baskerville (3) to predict the 
performance of a surface with a given geometry and to study the 

effect of varying geometry on the performance of enhanced surfaces. 
There are six basic geometric parameters for a triangular pitch 
circular finned tube bank: tube diameter; tube transverse and 



diagonal 


pitch 


; fin height, thickness 


, and 


spaci ng . 






It 


was 


found that for the 


same 


heat transfer (NTU) 


and 


pumping 


power, the volume of a 


heat 


exchanger 


decreases 


for 


smal ler 


tube 


diameter and tube 


pi tch . 


Opti mum 


values of 


f i n 


spaci ng 


and 


fin height were found 


for several 


specific 


f i n 
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thicknesses, tube diameters, and pitch. Thinner fins were found to 
be more effective in enhancing a surface than thicker fins for the 
same fin spacing in terms of fin thicknesses. 

This study extrapolated correlations based on experimental 
data for heat transfer and pressure drop beyond the range for 
which data exists. Experimental or theoretical work remains to be 
done to either prove the validity or correct any correlations that 

might be used in a similar parametric comparison study for the 
range of the geometric parameters where optimum performance 

appears to occur. 

Future studies on the geometry effects of heat exchangers in 
crossflow might be extended to geometries not limited to triangular 
pitch tubes, circular fins, or circular tubes. With the proper corre- 
lations parametric studies, such as this one, might be performed for 
all the other geometries a crossflow heat exchanger might assume. 
Use of a high-speed computer would allow plotting families of 

performance curves versus the various geometric parameters. 

Finally, performance might be defined in terms other than 

minimum volume for the same number of transfer units (NTU) and 
pumping power. A promising candidate might be minimum exchanger 
weight for the same NTU and pumping power. This would involve 
minimization of the volume of the material (tubes and fins) of an 
exchanger instead of the actual volume occupied by the exchanger 
for a given performance. Another dimension to this problem would 
be the materials from which the exchanger is constructed from. For 
example, although aluminum is lighter, copper has a higher heat 
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conductivity. However, materials are often dictated by service 
environment. Therefore, optimum performance for minimum volume of 
exchanger material would be a most useful measure to minimize the 
weight of an exchanger. This would insure the most efficient use of 
heat exchanger material on a per weight basis for a given 
performance. 
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APPENDIX I 

The purpose of this appendix is to present the derivation of 
several geometric relationships. 

A . Derivation of total heat transfer area on one side, A^, for 
circular finned tubes with fins of rectangular profile. 

Figure A1 shows the layout of a circular finned tube with fins 
of rectangular profile. 



The total heat transfer area including fin surface and edges 
of fins is: 



> 

H 

II 




[2(^d/a -,d// 4) + -,D^t,Oj/((SBp+.,)Dj; 


] (A1 ) 




Since, 








D = D 
f o 


4 2H,D^ 


(A2) 


and, 


2 2 

°,-°o 


= (D,-D )(D,+D ) = 4H,^D ^ + 4H,D ^ 
f o f o f o f o 


(A3) 


> 

H 

II 


[ 


1 + + ')/(Sgp + t,,)]] 


(A4) 


Letting : 










' H,/(Sgp - t,) 


(A5) 




ttLD [1 - 
o'- 


t- 2R^ (H^ + t^ + 1 ) ] 


(A6) 


and 


letting : 










= [1 + 2R^ (H^ + t^ + D] 


(A7) 


^ = 


.LDoGf 




(AS) 



B . Derivation of minimum flow area for gas, A^ . 

Figure A2 shows the triangular layout of a tube bank in 



crossf low . 



95 



FIGURE A1 

Circular finned tube with fins of rectangular profile. 
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FIGURE A2 

Triangular layout of a tube bank. 




96 



Minimum flow in the transverse space occurs when 
while minimum flow occurs in the diagonal space when 2a^ < a^. 
Therefore, minimum flow in the transverse space implies < 2X^-1 , 
and minimum flow in the diagonal space implies Xp<(X^+l)/2. In 
either case, Figure A3 shows fins reduce the minimum flow area. 

Therefore, in the calculation of minimum flow area, the fin 
cross sectional area must be subtracted. 

For Case I: Minimum free flow in the transverse space: 



and, 

A = L[X^.D - D - 2tp H.D /(S„^+tjD ] 
c To o f o f o BF f o 

= LD (Xt- - 2R,t, - 1) (A9) 

o T f f 

For Case I I : Minimum free flow in the diagonal space: 

Xj^ <(X^+1)/2 

and, 

A = L[2X_D - 2D - 4t,D /(S„^+tjD ] 
c Do o fofo BFfo 

= 2LD^[Xp - 2R^t^ - 1] (A10) 



C. Derivation of hydraulic diameter, D^^. 



For Case I: Minimum free flow in the transverse space: 



and, from equations (la), (A6), and (A9) 

^ 4(X^-2R^t^-1) X|_D^ 4X|_D^(X^-2R^t^-1 ) (All) 

^ " 7r[1+2R^(H^+t^+1 ) J " 

For Case I I : Minimum free flow in the diagonal space. 



Xp < (X^+1 )/2 
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FIGURE A3 

Layout of a tube bank with gas flow normal to paper. 
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and, from equations (la), (A6), and (A10) 

8(X -2R t.-l )X, D 8X, D (X_-2R,t,-1 ) (A12) 

U. _ U ff Lo _ Lou ff 

7r[1+2R^(H^+t^+1)J ■ 



D. Derivation of ratio of total heat transfer area on one side of the 



exchanger to total volume of the exchanger, /? . 



From equation (A6) and Figure 3: 

A-p 7 t[ 1 +2R, ( H^+t^+1 ) ] irGr 

R = — 1 = 1 — I — I = 1 — 

^ X^X, D X^X, D 

T T L o T L o 



(A13) 



E. 



Derivation of ratio of fin area to total transfer area, 






It is seen that; 



^ L[2[ 7t(D^V4) - 7t(D^^/4)] + 



f = 






= 27tLD R^(H. + t, + 1 + tV2H.) 

O T T T T T 



(A14) 



Therefore, 

A^ 2R^(H^ + + 1 + t^/2H^) 

A^ " LI + 2R^(H^t^+1)J 

= 2R^(H^ + + 1 + t^/2H^)/G^ (A15) 



F . Derivation of the angle define by the longitudinal centerline and 
the diagonal centerline, a ; see Figure 4. 

(A16) 



2 arctan 



X^-D /2 
T o 

X, D 

L L o J 



2 arctan (X^/2Xj^) 
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G . Derivation of the equivalent diameter used by Jameson (5) and 
Mirkovic (7) , 

The projected perimeter of a finned tube, is 2tt times the 

distance traveled axially along the tube and over the surface of 
the fins from one end of the tube to the other. (See Figure A1.) 
Therefore, 



2 7T L(2Hp + t^D + ) 

f o f o BF o 

f o BF o 



= 27tL(2R^ + 1) 



(A17) 



The equivalent diameter of a finned tube, d^, used by 
Jameson (5) and Mirkovic (7) is defined as twice the total air heat 
transfer surface area divided by tt times the projected perimeter of 
a tube. Therefore, 



2A^ D^[l+2R^(H^+yl)] 



TT I , 



2R^+1 



D 

o F 
2R +1 



(A18) 
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APPENDIX II 



The purpose of this appendix is to compare the method 
developed by Biery (4) as presented in Chapter I I I for predicting 
the Colburn j-modulus for gas flow normal to triangular pitch 
finned tube banks with the data given by Kays and London (2) and 
with several other correlations. 

Table A1 presented the percent deviation of the Colburn j- 
modulus predicted by Biery's method from the experimental value 
given by Kays and London (2) for the correspondi ng surface over 
the given range of Reynolds numbers, where 



Percent Deviation 



•*8 •*K&L X 1 00% 

•^K&L 



(A19) 



Table A1 shows the Colburn j-modulus predicted by Biery's 
method deviates from the experimental values by approximately 40 
percent for surfaces 93 and 94. However, Kays and London notes 

experimental uncertainty for heat transfer for surface 93 due to the 

necessity of estimating a contact resistance in the bi-metal tubes 
which were used. Another source of error for surfaces 92, 93 and 94 
is, according to the values given, Kays and London calculated the 
hydraulic diameter without accounting for the fins as outlined in 
Appendix I, Section C. However, the values given for hydraulic 

diameter for the other surfaces, which are data taken from Jameson 
(5), suggest the fins were accounted for in calculating the 

hydraulic diameter. Therefore, due to these errors, surfaces 93 and 
94 are discounted. The remaining body of data shows that Biery's 
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TABLE A1 

Deviation of the Biery (4) method for Colburn j-modulus in 
predicting experimental ly obtained value as presented by Kays and 
London (2), 



Surface"^ 


RE=1000 


RE=4000 


RE=7000 


RE=10000 


92 


-27.1% 


-18.5% 


-14.7% 


-12.2% 


93 


-41 .7% 


-40.7% 


-40.3% 


-40.0% 


94 


-37.3% 


-37.1% 


-37.0% 


-36.9% 


96 


-12.0% 


-18.8% 


-21 .4% 


-23.1% 


97A 


-16.3% 


-20.8% 


-22.6% 


-23.7% 


97 B 


-11.5% 


-15.5% 


-17.1% 


-18.0% 


98A 


+6.7% 


-3.6% 


-7.4% 


-9.7% 


98 B 


+17.8% 


+6.0% 


+ 1 .6% 


-1.1% 


98C 


+27.4% 


+13.0% 


+7.7% 


+4.5% 


98 D 


+6.3% 


-4.0% 


-7.9% 


-10.2% 


98 E 


+12.0% 


+ 1 .4% 


-2.5% 


-5.0% 


99A 


-10.1% 


-13.4% 


-14.7% 


-15.6% 


99B 


-1 .5% 


-4.3% 


-5.5% 


-6.2% 


*Surface 


designation 


corresponds 


to the 


Kays and London 



figure number. 
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method of predicting the Colburn j-modulus agrees with the 
experimental values given by Kays and London (2) to within an 
average per surface of 20% from the mean value with a mean 7% 
below the Kays and London data and with a standard deviation of 
8 %. 

The method given by Biery (4) for predicting Colburn 
j-modulus was also compared to correlations given by Briggs and 
Young (10) and Mirkovic (7). Caution should be exercised when 
comparing these correlations as the Reynolds number for each author 
Is based on a different diameter. To compare these correlations a 
factor must be included to convert the correlations to a common 
diameter. The correlation given by Briggs and Young is based on 
the outside tube diameter while the correlation given by Mirkovic is 
based on an equivalent diameter, d^, as presented in equation 
(A18). 



The 


Reynolds 


number used by Biery (4) 


is 


based 


on the 


hydrau 1 ic 


diameter 


given by Kays and London 


(2) 


and 


equat Ion 



(la). 

The correlation given by Mirkovic (7) when converted to a 
hydraulic diameter basis may be expressed as follows: 



•*M 



^Sl'^Pr 



2/3 



r.. 1.338 



224 


'l . 1 547X, D-D 
L o o 


.1 


1 

X 

r 

0 
o 

1 

o 

o 

1 


-.15 


'^bf' 


.25 


^h‘ 




D 

o 




D 

o 








d 

e 



RE-.338 

(A20) 



It should be noted that 
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S 



BF 




(A21) 



info 



which is the way Mirkovic originally included this factor in his 

correlation. It appears Mirkovic intended to develop a correlation 

which included fin thickness as a variable in this factor, although 

it is not. However, the Mirkovic correlation is slightly dependent 

upon fin thickness when expressed in terms of a hydraulic diameter 

due to the (D. /d ) term, 
h e 

Biery (4) noted the Mirkovic correlation gave very fine results 
when used to analyze the Mirkovic (7), Kays and London (2), and 
Jameson (5) data. However, this should be expected since the data 

produced the parameters and coefficients for the correlation. Biery 
reported the overall result in comparing the Mirkovic correlation 
with these data gave an average error of +.3 percent with a 

standard deviation of 3.6 percent. However, since the Mirkovic 

correlation is not directly a function of fin thickness, its more 

general validity outside the narrow range of data from which it 
was drawn is to be suspected. Also, the Mirkovic correlation 
includes a term that subtracts the diameter of the tube from the 
longitudinal distance between tube centers. In a triangular pitch 
tube bank, where the transverse tube pitch is large and the 
longitudinal tube pitch is small, the longitudinal distance between 
the tube centers can approach and sometimes be smaller than the 
diameter of the tube. 

The correlation given by Briggs and Young when converted to 



a hydraulic diameter basis may be expressed as follows: 
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•^B&Y “^ST^Pr ^ 





.1134 




.2 




^BF 




^BF 




_ ‘f _ 




_^f J 
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o 



.319 oe--319 
K t. 



(A22) 



Gianolio and Cut! (11) confirmed in a set of experiments that 
the Briggs and Young correlation predicted the air heat transfer 
coefficients obtained from their experiementa! data for induced draft 
for six tube-row banks, but had to include correction factors given 
by Ward and Young (12) for tube banks having a number of rows 

less than six. 

Biery*s method for predicting the Colburn j-modulus agreed 
with the Mirkovic and Briggs and Young correlations to within plus 
or minus 15 percent over the range of values of tube diameter, tube 
pitch, fin thickness, fin height, and fin spacing from which these 
correlations were based. However, these correlations deviated from 
Biery’s method by as much as 40 percent at extreme values outside 
the range of data. All three correlations were also in agreement by 

predicting an increase in the Colburn j-modulus with an increase in 

fin spacing or tube pitch or with a decrease in fin height or tube 

diameter. 

However, as fin thickness was increased, Biery's method and 
the Mirkovic correlation predicted an increase in the Colburn 
j-modulus, while the Briggs and Young correlation predicted a 
decrease. The experimental data taken by Mirkovic (7) presents an 
opportunity to examine variation of the Colburn j-modulus with fin 
thickness. In Mirkovic's run 12 and 13 all parameters were 
identical except for fin thickness. Mirkovic’s Figure 2 shows for 
run 13 (the thicker fins), the Colburn j-modulus is higher than for 
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run 12. Therefore, an increase in the Colburn j-modulus with 
increase in fin thickness seems to be the correct indication. 

In conclusion, all three correlations yield general agreement, 
except for variation of fin thickness. For this parameter two of the 
three correlations (Biery and Mirkovic) are in agreement both with 
each other and with a particular example from the Mirkovic data. 
Therefore, the choice of a correlation for predicting the Colburn 
j-modulus of a given surface, rests between the Biery method and 
the Mirkovic correlation. In the final analysis, Biery’s method was 
selected over the Mirkovic correlation. Biery’s method is more 
general and based on a larger body of data encompassing a larger 
range of geometries. Further, it is not restricted to geometries in 
which the longitudinal distance between tube centers must be 
greater than one tube diameter as in the Mirkovic correlation. 
Therefore, the Biery method for predicting the Colburn j-modulus is 
used later for predicting the performance parameters of a given 



surface. 
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APPENDIX I I I 

The purpose of this Appendix is to compare Jameson's (5) cor- 
relation for pressure drop across triangular pitch finned tube banks 
with the data given by Kays and London (2), 

As noted in Chapter IV, to compare the friction factor used in 
the Jameson correlation, fj, with the friction factor given in the 
data by Kays and London, a geometry factor must be included and 
the Reynolds number converted to a common diameter. With these 
changes the friction factor as defined by Kays and London (2) per 
equation (30) and as predicted by the Jameson correlation becomes: 



f = 1.532 







'Dj‘ 


X, D 






L o 




h 



-.25 



RE 



-.25 



(A23) 



Table All presents the percent deviation of the friction factor 
as predicted by the Jameson correlation, (f), and the experimental 
value given by Kays and London (2) for the corresponding surface 
over the given range of Reynolds numbers, where 

Percent Deviation = ^ 100% (A24) 

Table All shows the friction factor predicted by the Jameson 
correlation (f*) for surface 98D deviates from the experimental value 
by approximately 67%. Surface 98D is the one case of the surfaces 
listed in Table All which the minimum free flow area occurs in the 
diagonal spaces. Jameson (5) noted that his correlation did not 
accurately predict the friction factor for finned tubes in which the 
minimum free flow area occurs in the diagonal space. This 
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TABLE Al I 

Deviation of the Jameson (5) correlation for friction factor in 
predicting experimentally obtained value as presented in Kays and 
London (2) . 



Surface* 


RE=1000 


RE =4000 


RE =7000 


RE=10000 


92 


+33.5% 


+23.8% 


+20.0% 


+17.7% 


93 


+17.3% 


+9.7% 


+6.7% 


+4.9% 


94 


+9.9% 


+5.3% 


+3.5% 


+2.4% 


96 


+8.9% 


+9.8% 


+10.1% 


+10.4% 


97A 


+0.5% 


-0.1% 


-0.4% 


-0.5% 


97B 


-4.2% 


-1 .8% 


-0.8% 


-0.2% 


98A 


-2.1% 


+0.7% 


+ 1 .9% 


+2.7% 


98B 


-0.2% 


+0.3% 


+0.4% 


+0.5% 


98C 


+10.8% 


+10.3% 


+10.0% 


+9.9% 


98D 


+65.1% 


+67 . 2% 


+68.1% 


+68.6% 


98E 


+6.2% 


+7.8% 


+8.4% 


+8.9% 


99A 


+ 18.1% 


+17.7% 


+17.6% 


+17.5% 


99B 


+9.7% 


+10.6% 


+1 1 .0% 


+11.3% 



^Surface designation corresponds to the Kays and London (2) figure 



number . 
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conclusion appears to be supported by the results in Table All for 
surface 98D . Discounting this one surface, Table All shows the 
Jameson (5) correlation predicts the Kays and London (2) data to 
within approximately 20% of the mean value with a mean 7.8% above 
the Kays and London data and with a standard deviation of 7.3%. 
The Jameson correlation is used in Chapters V and VI to study the 
geometry effects of tube banks in which the minimum free flow 
area occurs in the spaces transverse to the air flow. 
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APPENDIX IV 



The purpose of this Appendix is to put together the equations 

and correlations introduced or derived in other Chapters and 

Appendices in an algorithm which may be programmed into a 

computer or calculator for predicting the performance parameters of 

a triangular pitch circular finned tube bank in crossflow given the 

geometry and the air and fin material properties. Hence, the only 

inputs required are tube diameter, tube pitch, fin thickness, fin 

height, fin spacing, and the air and fin material properties. 

The air and fin material properties are entered as a 

’’property constant” (P ) in the computation of fin efficiency as 

c 

outlined in Figure 5. 



Since: 




1.5 



( 10 ) 



m 



and : 




(A25) 



where from (3a): 



RE /X 



(A26) 



G 



c 



It can be shown: 



ml 



2 P j RE D 
c o 




(A27) 



where ; 



P 




(A28) 



no 



For dry air at 90^ F at atmospheric pressure = 1.17 10 ^ 
for aluminum fins and = 6.23 lO”^ for copper fins. 

The algorithm for predicting the performance parameters of a 
surface given the geometry and property constant is introduced in 
Figure A4. This algorithm calculates the performance parameters for 
a triangular pitch smooth tube and finned tube with radial fins of 
rectangular profile only. An algorithm that would include tapered 
fins, continuous fins, in-line tubes, and other various geometries 
could be written employing the methods outlined in this paper, but 
would require the applicable correlations and geometric relation- 
ships. It should be recalled that this algorithm is not accurate for 
predicting enhanced surface performance parameters for minimum 
free flow area occuring in the diagonal due to the inaccuracy of 
the Jameson (5) correlation for this case. However, the program 
should produce accurate results for both cases of minimum free flow 
for bare tubes. Therefore, the effect of varying tube pitch can be 
studied even for the case of minimum free flow occuring in the 
diagonal by using the performance parameters produced for smooth 
tubes. The program is written to produce performance parameters for 
bare tube Reynolds numbers 10000, 70000, and 130000. 

Values for only two Reynolds numbers are required to plot the 
performance parameters for a surface since the plot is a straight 
line on a log-log scale. However, a third value insures plotting 
accuracy. The range for REj^ of 10000 to 130000 was chosen since 
the typical value of RE^^/RE is 7 to 20. Therefore, performance 
parameters are produced for Reynolds number based on the 
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hydraulic diameter in the range of 500 to 18000. This is the 

approximate range most finned tube heat exchangers in crossflow 
operate and is the same range for which the Kays and London (2) 

data is given. 

The algorithm of Figure A4 was keyed into a T I -59 calculator. 
This calculator program is listed in Figure A5. Due to the length of 
the algorithm, the calculator program is divided into two halves. 
The first half calculates the parameters for the equations in the 

algorithm up to label A (Steps 000 to 422), and stores these values 
in data registers for use by the second half of the program. After 
the first half of the program has been run, the second half of the 
program (Steps 000 to 632) can be read in by magnetic card and 

run to generate the performance parameters. Table AIM lists the 

data registers in which the input and calculated parameters are 
stored. Where more than one parameter is listed indicates the 

corresponding data register is used for temporary storage of the 
indicated parameter while the program is running. 

Figure A6 shows a sample output generated by running the 

program of Figure A5. The program labels refer to the indicated 

parameter. The surface is defined in the output by a listing of the 
> 

input parameters. It should be noted that the only input parameter 
not listed in the output is P^. In the parametric study of Chapter 
VI the fin material was assumed to be copper and the fluid on the 
enhanced surface side was assumed to be dry air at 90° F for all 

cases. The property constant (P ) for these assumptions is 6.23 

10 '^. 
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The output of Figure A6 can be used to plot the performance 
parameters for the indicated surface. A comparison can then be 
made on the surface effectiveness as outlined in Chapter 2. Each 



surface 


generated in 


the 


parametric 


study of 


Chapter VI was 


compared 


on the basis 


of 


volume one 


side for 


the same pumping 



power and number of transfer units (Case C comparison of Chapter 

II). Each surface was compared to a bare tube surface where X-^. = 

Xp = 2 and = 1”. To standardize this comparison the 45 degree 

line (slope = 1 on log-log scale) was drawn from the point on the 

performance parameter plot of the reference bare tube surface where 

log (NTU/V ) = 2.0. 
os 

The 45 degree line described above and the plot of the 
performance parameters of a surface are straight lines. Therefore, 
it should be finally mentioned that as an alternative to plotting a 
performance parameter to make a volume comparison to another 
surface, linear equations may be written and solved simultaneously. 
This permits an analytical comparison which is faster if in the form 
of a program and also more accurate. 
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INPUT 

^T’ ^D’ ^BF’ 

P , D (inches) 
c o 



Calculation of Various Parameters for Later Use: 

X|_ - (Xq^ - 

^ " *f> 

= 1 + 2R^(t^ + + 1 ) 

^ 

A^/A^ = 2R^(H^ + + 1 + t^/2H^)/G^ 

NrI " 

^HD = - 2Rftf - DX^/- 

= 16/ 

RE|^/RE = ^X^X|_D^Tr 

d = D G,/(2R, + 1 ) 

e off 



°J - 



(H/2Sbf) 



i(X D /d -1) 2 + 1(X D /d -1)~2 
I o e Doe 



In (r /r ) = - In (1 + 2H-) 
o e T 

C. = -.00259 - .0084 In (r /r ) 

^ o e 

D. = -.02426 - .4521 In (r /r ) 

^ o e 



Figure A4: Algorithm for Prediction of Performance Parameters for 
Triangular Pitch Circular Finned Tube Surfaces and Bare 
Tube Surfaces. 
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Figure A4. (Cent) 
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Figure A5. T I -59 Calculator Program for Predicting Performance 
Parameters (First Half). 
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520 55 - 


560 95 = 


6 U U 43 R G L 


521 43 RCL 


561 65 X 


601 00 00 


d22 2U 2u 


562 43 RCL 


602 55 


523 95 = 


.563 14 14 


603 43 RCL 


524 42 STD 


564 " 65-..X 


604 31 31 


525 36 36 


565 43 RCL 


605 54 > 


526 23 LNX 


0 b b b b b' b 


6 0 6 45 Y X 


D 2 6 0 


56? 55 ^ 


607 93 


528 43 RCL 


568 43 RCL 


608 04 4 


529 29 29 


569 31 31 


6= U 9 9 4 ~ 


530 94 +•••■- 


570 33 X£ 


6 1 0 95 = 


531 85 + 


571 95 = 


611 65 X 


532 43 RCL 


572 28 LOG 


612 43 RCL 


533 28 28 


573 69 DP 


613 14 14 


534 95 = 


•5 r " 4 0 6 0 6 


614 65 X 


535 94 


— — — 

• J • J - Z: 


615 43 kGL 


536 85 + 


576 03 3 


Q 1 "O 3 O -I' 


537 43 RCL 


577 05 5 


617 55 ^ 


538 36 36 


5 7 Q 6 5 X 


618 43 RCL 


539 95 = 


579 43 RCL 


619 31 31 


540 94 +/- 


580 14 14 


620 33 K- 


541 85 + 


581 45 YX 


621 95 = 


542 01 1 


-;l:v Q--; 


622 28 LOG 


543 95 = 


5 8 3 03 3 


623 99 PRT 


544 65 K 


584 09 9 


6 3 4 U 6 6 


545 43 RCL 


585 06 6 


625 0 0 0 


546 18 18 


b y b y 4 ^ ~ 


626 00 0 
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TABLE AMI 

Data Registers used for Storage of the Indicated Parameters 
for the Program in Figure A5. 



Data Register 


Parameter (s) 


Data Register 


Parameter ( s ) 


00 


D 

o 


23 


d 

e 


06 


P 

c 


24 


"j 


07 




25 


E 


08 


^BF 


26 


"e 


09 


‘f 


27 




10 


H. 

T 


28 


Ca 


1 1 


Xj 


29 


Da 


12 


Xl 


30 


°h 


13 


Xd 


31 




14 




32 


N e 


15 




33 


A^-b/V 
F os 


16 




34 


ml _ 

^ ’ f-Corr 


17 


/? 


35 


C C e'^‘ 
^T’ D’ 


18 


Af/Ay 




log(NTU^/V^^) 


19 


^HT’ -^'e 


36 


's> ^-L 


20 


^ Xj ^4 ’ ^ ’ 


37 


RE 




log(NTU /V ) 

^ S OS 


38 


log(P^/V^s> 


21 


N,,Ch 


39 


'°9(Ps/^os^ 


22 


RE^/RE 
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1 . 

0 = 02 
0 , 02 
0 = 5 



1 U000 = 




!-'0000 = 
13 = 59794349 
12 = 12723271 
3 = 287570581 
2 = 043031614 



14 = 3372681 
12 = 87685563 
3 = 437757844 
2 = 205414183 




R 

P 



RE 




N 

P /V @RE,=10000 
e os N 

P /V @RE.. =10000 
s os N 

NTU /V @RE. =10000 
e os N 

NTU /V @RE^, = 10000 
s os N 

etc . 



Figure A6 . Sample Output 



Generated by Program listed in Figure A5. 



1S9468 



Thesis 

R6429 Robinson 

Geometry effects on 
triangular pitch cir- 
cular finned heat 
exchangers in cross- 
flow. 




Thesis 

R6429 

c,l 



199468 

Robinson 

Geometry effects on 
triangular pitch cir- 
cular finned heat 
exchangers in cross- 

flow. 



